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ABSTRACT 

The transverse motions of nearby dwarf spheroidal (dSph) galaxies contribute line-of-sight com- 
ponents that increase with angular distance from the dSph centers, inducing detectable gradients 
in stellar redshift. In the absence of an intrinsic velocity gradient (e.g., due to rotation or stream- 
ing), an observed gradient in the heliocentric rest frame (HRF) relates simply to a dSph's systemic 
proper motion (PM). Kinematic samples for the Milky Way's brightest dSph satellites are now suffi- 
ciently large that we can use stellar redshifts to constrain systemic PMs independently of astrometric 
data. Data from our Michigan/MIKE Fiber System (MMFS) Survey reveal significant HRF veloc- 
ity gradients in Carina, Fornax and Sculptor, and no significant gradient in Sextans. Assuming 
there arc no intrinsic gradients, the data provide a relatively tight constraint on the PM of Fornax, 
ifJ-a^^^fJ-s^'^) = (+48 ± 15, -25 ± 14) mas century"^ that agrees with pubhshed HST astrometric 
measurements. Smaller data sets yield weaker constraints in the remaining galaxies, but our Carina 
measurement, (^^ A^a^^^) = (+25 ±36, 4-16 ±43) mas century"^, agrees with the published astro- 
metric value. The disagreement of our Sculptor measurement, [iJ-a^^ t i^f^^) = (—40 ±29, —69 ±47) 
mas centur y"^, with astrometric m easurements is expected if Sculptor has a rotational component as 
reported bv lBattaglia et all (120081 ). For Sextans, which at present lacks an astrometric measurement, 
we measure (/^^ IJ-^^^) = (—26 ± 41, +10 ± 44) mas century"^. 

Subject headings: galaxies: dwarf — (galaxies:) Local Group — galaxies: individual (Carina, Fornax, 
Sculptor, Sextans) 



1. INTRODUCTION 

In recent years there has been an explosion in the 
quantity of kinematic data available for individual stars 
in the brightest Local Group dwarf spheroidal (dSph) 
; galaxies (e.g., Battag ha et al.l 120061: i Koch et al.l 120071 : 
• iMateo et al1l2008i: ,Walker et al.ll2008aD . Many of the ma- 
jor surveys producing such data are designed to quantify 
the amount and distribution of dark matter in these sys- 
tems, but the improved statistics benefit other investiga- 
tions as well. We now know that what were once thought 
to be simple dSph stellar populations often contain 
multiple compone nts separable by age, metallicity an d 
■ kinematics (Tolst oy et al.l 120041: IBattaglia eFaT '2008D. 
Some dSphs display evidence for lo calized kinematic 
substructure (jKlevna et a"ni2003t IWalker et all i2006bl ). 
while some exhibi t tidal streaming motion among their 
outermost stars jMuiioz et all 120061 : iSohn et al.ri2007l : 
IMateo et al.l 120081 ). Taken together, all these results in- 
form models of galaxy evolution at the smallest scales. 

Here we put the available kinematic data sets to 
another new use — measuring systemic proper motions 
(PMs) of dSphs ind e pende ntly of astrometric data. 
iKaplinghat fc Strigaril (|200^ have recently argued that 
one can use line-of-sight (LOS) velocity samples of > 
1000 stars to measure PM with precision similar to that 
which is achieved with HST astrometry, while samples of 
> 5000 can do significantly better. The technique relies 
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on the detection of "perspective rotation" : at large angu- 
lar radii, a dSph's transverse motion has a non-negligible 
component along the LOS, causing spectra to become in- 
creasingly redshifted along the direction of the PM vec- 
tor. In the absence of an intrinsic velocity gradient (e.g., 
due to true rotation and/or tidal streaming), the strength 
and orientation of an observed gradient relates simply 
to the systemic PM. This effect is frequently consid- 
ered in kinematic studies of the Magellanic Clouds (e.g., 
iFeast et al.lllQeiHvan der Marel et al.ll2002f ). and has re- 
cently been used to cons train the PM of M31 from the ve- 
locities of its satellites ()van der Marel fc Guhathakurtal 
I2OOI . 

Here we apply the perspective technique for the first 
time to large dSph velocity samples and present the 
resulting constraints on the PMs of Carina, Fornax, 
Sculptor and Sextans. The data result from our Mag- 
ellan/MMFS observ ations as of 2008 A ugust and are 
presented elsewhere (jWalker et al.|[2008al Paper II). To 
our MMFS d ata we add the Fornax kinematic data of 
IWalker et al.l (l2006a), which contribute 155 stars not ob- 
served with MMFS. We do not combine with other pub- 
lished data, but have verified that, with one exception 
that owes to different spatial sampling (see Section [5] 
discussion of Carina), doing so would leave our results 
qualitatively unchanged. The samples contain velocity 
measurements for the following numbers of stars: Ca- 
rina: 1982 stars (774 members); Fornax 2793 stars (2610 
members); Sculptor: 1541 stars (1365 members); Sex- 
tans: 947 stars (441 members). 

2. VELOCITY GRADIENTS 

We first demonstrate that the available kinematic data 
exhibit significant velocity gradients. For each dSph, Pa- 
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per II presents LOS stellar velocities, V, measured in the 
heliocentric rest frame (HRF). Allowing for an HRF ve- 
locity gradient k = dV/dR', where R' is angular distance 
from the dSph center in the direction of the gradient, the 
data have likelihood 
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where is the measurement error, cvb is the internal 
velocity dispersion, and 9i, 6o are position angles of the 
star and gradient, respectively. Walk er et alJ ()2008bL Pa- 
per III hereafter) introduce an algorithm that evaluates 
the membership probability, Pm, of each star according 
to its velocity, magnesium index and position. Assigning 
weights to the data points according to their membership 
probabilities, the expected log-likelihood is 
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For a given 9q, the estimates'' {V), ayg and k t ake the 
values that maximize i5(lnL). Following Walke r^t al.l 
(|2008b), we obtain estimates by setting equal to zero 
the partial derivative of E{hiL) with respect to each pa- 
rameter, and then solving iteratively. For example, k is 
calculated in each iteration as 
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For each dSph we consider possible position angles for 
the velocity gradient of 6*0 = {0°, 3°, 6°, 180°}. We 

quantify the significance of the maximum gradient, k^ax, 
via Monte Carlo simulation using permutations of the 
real {Vi, crv. , i?^, 9i) data. In each of 1000 realizations we 
re-assign [Vi, aVi, PmA tuples randomly to 9i) pairs. 
Thus the simulated data sets have the same spatial sam- 
pling and overall velocity distribution as the real data, 
but scramble any existing correlation between velocity 
and position. We define the significance, pj, , of the 
maximum velocity gradient measured from the real data 
to be the fraction of simulated data sets that fail to pro- 
duce, at any position angle, a gradient as large as the 
maximum gradient observed in the real data. 

Table[l]lists for each dSph the values of kmax and pf, 
corresponding to the position angle, 6'o,„„^ , that gives the 
largest velocity gradient. Carina, Fornax and Sculptor all 
exhibit HRF velocity gradients that are significant at the 
Pf, > 0.973 level. Sextans shows no significant gradi- 
ent, with p^ — 0.753. Despite their significance, none 
of the gradients inflate estimates of the global velocity 

* X denotes the estimate of X. 



dispersion — estimates (V) and avo are identical to those 
obtained in the analysis of Paper III, where we assume 
no velocity gradient. 

3. DSPH PROPER MOTIONS 

The observed velocity gradients can result from the 
perspective effect and/or some combination of intrinsic 
rotation and streaming motions. Here we assume that 
the latter two phenomena are negligible, in which case a 
gradient reflects the transverse motion of the dSph. Sup- 
pose Vrei{(x,S) is the projection of the relative motion 
between Sun and dSph along the line of sight defined by 
equatorial coordinates {a, 6). To an observer instanta- 
neously at the Sun but comoving with the dSph — i.e., in 
the dSph rest frame (DRF) — a dSph star at (a, S) has 
LOS velocity 



Vdrf = V - Vrel{a,S), 
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where V is the star's HRF velocity. In Appendix |X] we 
derive formulae that express Vrei{a,S) in terms of the 
components {iia, fJ-s) of the dSph's systemic PM. 

We assume a given dSph has a Gaussian velocity dis- 
tribution with mean (V") drf = and variance ay , with 
no intrinsic gradient. The HRF gradient (Section[2|) then 
reflects the smooth variation in Vrei across the face of the 
dSph, and the likelihood function becomes 
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Again assigning weights to each star according to prob- 
ability of dSph membership, we obtain the expected log- 
likelihood 
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We measure PM by maximizing E{h\L) over the pa- 
rameter pair {^JLg , ^s ) , wh ich specifies Urei(Q!i, <5i) accord- 
ing to Equations lA3l - [AEl During this procedure we hold 
the velocity dispersion trvb fixed at the value that maxi- 
mized the expected log-likelihood in Equation [31 These 
values are ctvo = 6.6 km s~^ (Carina), 11.6 km s~^ (For- 
nax), 9.2 km s~^ (Sculptor) and 7.9 km s~^ (Sextans). 

Figure [1] displays and Table [1] lists the resulting PM 
measurements. Contours in the fJ-a^^^l^f^^ plane of 
Figure [T] enclose regions containing 68%, 95% and 99%, 
respectively, of the volume underneath the exp[E{lnL)] 
surface. Points with errorbars indicate published, as- 
trometric PM me a surements (j Schwcitzer et al. 19951 
iPiatek et aL| 120031 : iDinescu et al . 2004; Piatck et al] 
11006, 200?^^^ 

For two (Carina and Fornax) of the three dSphs with 
astrometric PMs, the PMs derived from MMFS data 
agree well with astrometric values. We obtain our tight- 
est constraint on the PM of Fornax, for which we measure 
(Ai^^^,/if^^) = i+AS ± 15, -25 ± 14) mas century-^. 
The la error el lipse in Figure [T] ove rlaps the measure- 
ments of both IDinescu et al.l ()2004[ ) and IPiatek et al.l 
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(|2007f ). It should be noted that the Fornax measure- 
ment bv lPiatek et alj (|2007f ). based on HST astrometry, 
has the best precision of any dSph PM measurement (±5 
mas century"^). It is enco uraKing that our me asurement 
of Fornax agrees with the iPiatek et al.l (|2007f) value and 
has precision similar to that of HST measurements for 
other dSphs. 

Our measurement for Sculpt or disagrees with the 
two published ast rometric PMs (jSchweitzer et al.l 119951 : 
IPiatek et al.ll200^ . which also disagree with each other. 
However, if either of the astrometric PMs is correct, dis- 
agreement with our result is to be expected. If we use 
Equation |4] to remove the perspective effect due to ei- 
ther of the astrometric PMs, the resulting Sculptor ve- 
locities exhibit a strong residual DRF velocity gradient 
(fc > 6.0 km s-Meg-i; pj^^^^ > 0.994). Thus, if ei- 
ther of the astrometric PMs is correct, our velocity data 
imply that Sculptor has a strong intrinsic velocity gra- 
dient. Wheth er attributable to a ro tational component, 
as argued by I Battaglia et al.l (120081 ) on the basis of in- 
dependent kinematic data, or to tidal streaming, such 
a gradient would invalidate the assumptions inherent to 
our PM measurement. 

4. DISCUSSION 

We have presented the first constraints on dSph PMs 
that are derived from stellar redshifts. Given the sizes 
of our kinematic samples, we are able to measure PM 
with appr oximately the precision predi cted from simu- 
lations bv iKaplinghat fc Strigaril (|2008l K08 hereafter). 
K08 adopt a more sophisticated model for a dSph's in- 
trinsic velocity distribution, marginalizing over several 
parameters that describe the gravitational potential and 
orbital anisotropy. However, our simple assumption that 
the dSph velocities are drawn everywhere from a sin- 
gle, isotropic Gaussian distribution is consistent with 
the fiat em pirical velocity dispersion profiles of dSphs 
([Walker et al. 2007), and yields similar constraints as the 
K08 method. 

IBattaglia et al.l (|2008f ) conclude that Sculptor rotates, 
based on their detection of a significant LOS velocity gra- 
dient in the Milky Way rest frame. We detect a similar 
gradient, but note that it can be accounted for entirely 
by perspective effects, provided that Sculptor's PM dif- 
fers from both of the astrometric measurements. In order 
to examine the implications of our PM measurement for 
a non-rotating Sculptor's orbit, we have integrated the 
orbit backward in time for 6 Gyr using code provided by 
Slawomir Piatek. Assuming a static MW potential with 
disk, spheroid and halo compon ents identical to those 
adopted bv lJohnston et al.l (|l999[ ). our measurement im- 
plies a plausible Sculptor orbit {rperi = 66^45 kpc, Vapo — 
149^^03 kpc) with eccentricity similar to that implied by 
the HST astrometric PM {rpe„ = &8lLlT,rapo = 122125 
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Fig. 1.- — - Proper motions measured from stellar redshift data. Con- 
tours indicate Icr, 2cr and 3cr constraints on the PM of each dSph, ob- 
tained by evaluating Equation[7]using the MMFS velocity data. Points 
with error bars indicate pu blished PM measur ements de rived from pre- 
cise astrometry (Carina — IPiatek et al.| 120031 ; Fornax— jPinescu et all 
200l|Platek et al J 20071 : Sculptor— [Schweitzer et al.|1995l ; |Piatek et al.| 
200€ )■ Where two measurements exist, the HST measurement has 
smaller error bars. 

kpc; IPiatek et al1l2006[ ). 

Conclusions regarding Sculptor's status then depend 
on the initial assumption. If one assumes and corrects 
for (either of) the published astrometric PMs, the sig- 
nificant residual DRF velocity gradient implies a strong 
intrinsic gradient that is likely due to so me combination 
of ro tation and/or tidal streaming fsee IBattaglia et al.l 
12001 . If one assumes instead that Sculptor has no in- 
trinsic velocity gradient, then the HRF gradient implies 
a PM that is inconsistent with either of the published 
a strometric measurem ent s . 

iMunoz et al.l (|2006f l detect a strong LOS velocity gra- 
dient from Carina members projected at large angular 
distances {R > 1°), well outside the region sampled by 
MMFS. Had we combined their extended kinematic sam- 
ple with our MMFS data, we would have measured a 
PM of (/^^^^, ^^) = (+56 ± 29, -1 ± 33), in poorer 
agreem ent with astrometric measurements. IMufioz et al.l 
(|2006) argue convincingly that the outermost regions of 
Carina exhibit tidal streaming motions, which would in- 
validate our PM measurement that includes their outer 
sample. It is encouraging that in the inner regions where 
no evidence for tides exists, our PM measurement stands 
in excellent agreement with the astrometric measure- 
ment. In any case. Carina and Sculptor help illustrate 
that PMs measured from stellar redshifts ought to be 
regarded with caution, as they are only as valid as the 
assumption of no intrinsic gradient. The agreement, to 
high precision, of our Fornax PM with the astrometric 
measurement serves as strong evidence that Fornax has 
no intrinsic velocity gradient over the kinematically sam- 
pled region. 

Our results include the first PM measurement of any 



TABLE 1 

HRF Velocity Gradients and Proper Motions 
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i/s/deg) (dcg) (mas/cent) (mas/cent) 

2.5 ±0.8 21 0.973 +25 ± 36 +16 ±43 

3.3 + 0.2 120 > 0.999 +48 + 15 -25 + 14 

5.5 + 0.5 21 > 0.999 -40 + 29 -69 + 47 ^ -50 

2.1 + 0.8 120 0.753 -26 + 41 +10 + 44 
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kind that has been made for Sextans. Our measure- 
ment implies that Sextans, at a distance of 86 kpc 
(jMateol [19981 ) ■ is receding from its perigalacticon dis- 
tance of rperi — 66lg^ kpc toward an apogalactic dis- 
tance of Tapo = 129133^. We note, however, that the 
relatively large error bars on Sextans' PM accommodate 
orbital eccentricities ranging from 0.25 — 0.89 within the 
95% confidence interval. The most radial of these or- 
bits would bring Sextans within ~ 5 kpc of the Galac- 
tic center. If one takes the most likely PM at face 
value, it becomes unlikely that Sextans is a member 
of a stream associated with other well-known objects 
in the Milky Way Halo. None of the PMs predicted 



by iLvnden-Bell fc Lvnden-Belil (|1995f ) for possible Sex- 
tans associations with Sculptor, Sculptor/Fornax, or Pal 
3/Pal 2/NGC 5824 falls within even the large error bars 
of the measured PM (although none of the stream asso- 
ciations is ruled out beyond the 2a level). 
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MM acknowledges NSF grants 0206081, 0507453, and 
0808043. MOW acknowledges support from the STFC- 
funded Galaxy Formation and Evolution programme at 
the Institute of Astronomy, University of Cambridge. 



APPENDIX 
DWARF REST FRAME (DRF) 

Suppose Ajj and A* are the (three dimensional) position vectors of a dSph and one of its stars, respectively, specified 
in a coordinate system with origin at the Sun. The relative motion between Sun and dSph along the line of sight to a 
star whose equatorial coordinates are (a*, (5,) is the scalar projection 

A, 

Vrei{a^,S^) ^ ■ {Ad), (A1) 

where A = dA/dt. 

In order to evaluate Equation lAll we define Cartesian coordinates such that the -l-a;, —y, and +z axes point toward 
(a, 6) — {6'^ , 0'^'^^) , toward the vernal equinox, and toward the North Celestial Pole, respectively. Thus the position of 
a star is fully specified by A, = (^*x'^ + ^*yy + ^^'^ the unit vector B, = A*/A, has components 

-B*^ = cos((5») sin(a,); (A2) 
B^^ = - cos((5*) cos(a,); (A3) 
B,^ = sin(,5,)- (A4) 

The (three dimensional) velocity of a dSph at distance Ad , with HRF line-of-sight velocity Vd — Ad and HRF proper 
motion (naj^is) has components 

Ad^ — Vd cos{Sd) sin{aD) + Ad^J-o cos{Sd) cos(q!d) — ADfJ-s sui^Sd) sin(Q;£)); (A-5) 

ADy = -Vd cos((5d) cos(aD) -I- AD^J.s sin((5£,) cos(q;d) + ^dMq cos{6d) sin(a£,); (A6) 

A-D, = Vd sm{dD) + ADfJ^s cos{6d)- (A7) 
Thus the relative motion between dSph and Sun projects along the star's line of sight with magnitude 

Vrei ia^,6^) = B^ ■ Ad = B.,^Ad, + B^^ Ad^ + B.,^ Ad, ■ (A8) 
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